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> I 1 INTRODUCTION 



ABSTRACT 

We examine the properties of satellites found in high resolution simulations of the 
local group. We use constrained simulations designed to reproduce the main dynamical 
features that characterize the local neighborhood, i.e. within tens of Mpc around the 
Local Group (LG). Specifically, a LG-like object is found located within the 'correct' 
dynamical environment and consisting of three main objects which are associated with 
the Milky Way, M31 and M33. By running two simulations of this LG from identical 
initial conditions - one with and one without baryons modeled hydrodynamically - 
we can quantify the effect of gas physics on the z = population of subhaloes in 
an environment similar to our own. We find that above a certain mass cut, M su b > 
2 x 10 8 h~ 1 M Q subhaloes in hydrodynamic simulations are more radially concentrated 
than those in simulations with out gas. This is caused by the collapse of baryons 
into stars that typically sit in the central regions of subhaloes, making them denser. 
The increased central density of such a subhalo, results in less mass loss due to tidal 
stripping than the same subhalo simulated with only dark matter. The increased mass 
in hydrodynamic subhaloes with respect to dark matter ones, causes dynamical friction 
to be more effective, dragging the subhalo towards the centre of the host. This results 
in these subhaloes being effectively more radially concentrated then their dark matter 
counterparts. 



Bullock et al.ll2000l; ISpringel et al.ll200ll ; IStoehr et al]|2002l ; 
Gao et al.ll2004l ; iGill et alj|2004 ). 



In the current paradigm of structure formation, large ob- 
jects such as galaxies or clusters, are believed to form hier - 
archically, or through a "bottom- up" (1 White fc Rees|[l978l ) 
process of merging. Over dense peaks in the primordial 
density field collapse at high redshift due to gravitational 
instability and then merge with each other to form suc- 
cessively larger objects. N-body simulations of this pro- 
cess originally produced fairly smooth cold dark matter 
(CDM) halos jPrenk et al.lll985l ; iDavis et alJll985h homog- 
enized by poor resolution. More recently, studies exploit- 
ing high resolution simulations have managed to resolve a 
wealth of sub-galaxy sized substructures (subhaloes) that 
can survive the vio l ent processes assoc i ated with merging 
iGhigna et a.1.1 Il998l : iKTypin et al.lll999l ; iMoore et al.lll999l; 



High resolution numerical simulations that resolve sub- 
structures can be computationally expensive, and many au- 
thors choose to simulate just the gravitati onal forces, mod- 
eling the baryons semi-analytically (e.g. ICole et all |2000| ; 
iBower et alj|2006h . Alternatively, baryon s may be modeled 
hydr o dynamically dKatz fc Gunnl 1991 ; I Navarro fc Benz 
I l99l[ : iKatd 1 19921 ; INavarro fc White! Il994l ; iKravtsov et~ 
Il997f ). with the advantage that the thermal state of the 
intergalactic medium may be directly calculated. The in- 
clusion of baryons affe cts the background dark matter: 
iBlumenthal et ah! |l986) showed that dissipative baryons 
will lead directly to the ad iabatic contraction of the halo 
(see also iGnedin et al.|[20o3 ) and are thus a critical ingredi- 
ent in determining halo and subhalo properties. 
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More recently a number of authors have exam- 
ined these effects - gas dissipation a nd star formation 
- on galaxy and cluster s ize d hal os feailin et al.l 120051: 
Sommer-Larsen fc Limousin! |2009| ; | Romano-Diaz et al.l 



20091 ). A number of studies have also looked at the effect 
on substructure. iNaeai fc Kravtsovl (|2005h noted that 
the radial distribution of substructures in Adaptive Mesh 
Refinement simulations was less concentrated than both 
the radial di stribution of subhaloes in dark matt er only 
runs (see also lDiemand et~aL 2004 ; Gao et a l. 2004) as well 
as the dark matter itself. llNFagai fc Kravtsov 1 20051 ) ruled 
out numerical resolution as the cause showing that this bias 
was caused by the tidal stripping (of dark matter) from a 
subhalo's outer regions: the bias disappears if subhaloes are 
chosen according to their mass at accretion. Furthermore, 
IMaccio etT al. (2006) simulated a Galactic mass halo twice - 
once with baryons using Smooth Particle Hydrodynamics 
(SPH) and once with just dark matter in order to study 
the effect of baryons on lensing. They found that owing 
to adiabatic contraction in the dense cores of subhaloes, 
their SPH run produced an over abundance of subhaloes 
in the inner regions of the halo as well as an increase by 
a factor of 2 in the absolute number of subhaloes with 
respec t to the dark matter run. Finally, [Weinber g et al.l 
(2008) examined the effect of baryons on substructures 
of group and cluster sized (10 12 — 10 14 M©) halos with a 
more statistical approach. By running the same simulation 
with and without gas dynamics, they found an overall 
good agreement between the positions and mass function 
of subhaloes across a large dynamical range with the two 
different techniques. Like previous studies, they also found 
that the baryons steepen the potential well of both the 
host and subhaloes, making the core of the substructure 
more resistant to tidal stripping, resulting in the relative 
depletion of dark matter subhaloes towards the centre. 

Baryonic substructure may effect the central DM dis- 
tribution in a different way. Satellites crossing the cen- 
tral part of their host halo loose energy due to dynamical 
friction. This energy is being deposited in the DM com- 
ponent, leading to 'heating' of the DM and consequently 
this may lead to the soft ening of the central DM cusp 
||E1-Zant et alj l200ll , |2004| ). This has been supported by 



the fu ll galaxy formation sim u lation s of I Rom ano- Dia z et al 



2008); TRomano-Diaz et all l|2009h and Johansson et al 



2009) 



One of the main new feature of the present study is 
that it is based on constrained simulations of the local uni- 
verse, where 'local' means a few tens of Mpc around the LG. 
We designed our simulation's initial conditions by impos- 
ing constraints derived from observational data of the LG's 
environment; specifically we use surveys of galaxy peculiar 
velocities and the distribution of nearby X-ray selected clus- 
ters of galaxies. The result is a simulations of a 'LG' that sits 
within the proper dynamical environment. This includes the 



correct motion and position of the Local Supercluster, the 
Virgo and Coma clusters and the Great Attractor. The main 
virtue in using constrained simulations is that it provides a 
numerical laboratory for testing the formation of the local 
universe, and in particular the Local Group, in a numerical 
environment that closely matches our actual neighborhood. 
Thus, much of the cosmic scatter involved in comparing un- 
constrained simulations with specific cosmological objects, 
is removed. 

Other studies typically exa mine the effect of baryons 
in either cluster sized halos (|Nagai fc Kravtsovl 120051 ; 
IWeinberg et alj 120081 ) or i n a random Milky Way sized 



halo I Maccio et al.l 20061). We h a ve extended previ- 
ous work dNagai fc Kravtsovl 120051 ; IMaccio et ail I2OO6I ; 
IWeinberg et al-lboOSl ) by examining with improved mass and 
force resolution the the effect of baryons on substructures in 
three galactic halos that resemble the Local Group. 

We study how the effect of dynamical friction is affected 
by the relative difference in the tidal mass loss of a substruc- 
ture with and without gas-dynamics. By using two simula- 
tions - one with and one without gas dynamics - we are able 
to determine where and how the trajectories and properties 
of the same subhalo in the two different runs diverges. Our 
simulations also have unprecedented mass and force resolu- 
tion, an improvement on previous work of between 1 and 
2 orders of magnitude. We have tuned our simulations to 
make it very easy to identify matching subhaloes in both 
simulations. We thus follow in detail how the stripping of 
an individual subhalo is effected by the presence of baryons 
and how the stripped subhalo is then affected by dynamical 
friction. 

This paper is organized as follows. In Sec. [2] we describe 
our methods, specifically our simulations (Sec. l2.2| ). our halo 
and subhalo finding technique (Sec. I2.3|l and how we match 
dark matter subhaloes to SPH ones ("Sec. I2.4")l . We present 
our results in Sec. [3] and conclude in Sec. 3] 



2 METHODS 

In this section we describe how the simulation were con- 
strained and run, the method we have employed to identify 
substructures in each run and how we match SPH subhaloes 
to DM ones. 



2.1 Initial Conditions 

In this section we describe how we generated our constrained 
initial conditions. The simulations used in this paper are 
part of the CLUES (Constrained Local UniversE Simula- 
tions) projeclQ whose main goal is to run simulations that 
mimic as accurately as permitted by current observations, 



1 http://clues-project.org 
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Figure 1. The 2 = 0, DM (left), gas (centre) and stellar (right) distributions in our simulation of the Local group. We smooth the 
density field at each point and plot the Local Group in a box with side length is 1.3 h~ x Mpc. We identify as our "Milky Way" the galaxy 
at the top of the box, "M31" as the galaxy towards the centre and "M33" as the galaxy at bottom of the box. This box sits at roughly 
the correct distance from both Virgo and Coma mass clusters, as well as an object similar to the Great Attractor. 



the formation of the Local Universe from cosmological initial 
conditions. 

The following procedure is employed in generating con- 
strained simulations of the LG. We use as observational 
constraints peculiar veloci ties drawn from t he MARK III 
dWillick et al.lll997T). SBF ilTonrv et alj|200ll ) and the local 
volume l)Karachentsev et al J 12004 ) galaxy catalogs and the 
position and vir ial properties of nearby X-r ay selected clus- 
ters of galaxies llReiprich fc Bohriri gcr 2002). The Hoffman- 
Ribak algorithm (Hoffman fc Ribak)ll99lf T is used to gener- 
ate the initial conditions as constrained realizations of Gaus- 
sian random fields on a 256 3 mesh. These are used to per- 
form a few low resolution DM-only N-body simulations. It 
should be noted that the environmental restrictions do not 
constrain the properties (e.g. position relative to eachother, 
mass, etc) of LG members explicitly. Hence the low reso- 
lution simulations are searched and the one with the most 
suitable LG-like object is chosen for follow up, high reso- 
lution re-simulation. High resolution extension of the low 
resolution constrained realizations are obtained by creat- 
ing an unconstrained realization at the desired resolution, 
FFT-transforming it to k-space and substituting the uncon- 
strained low k modes with the constrained ones. The result- 
ing realization is made of unconstrained high k modes and 
constrained low k ones. 



2.2 The Simulations 

In this section we explain details of the simulation run and, 
for the gas dynamical simulation, the rules we employed for 
star formation and feedback. We choose to run our simula- 
tions using standard ACDM initial conditions, that assume 



a WMAP3 cosmology l|Spergel et alJl2007t >. i.e. fl m = 0.24, 
fif, = 0.042, S7a = 0.76. We use a normalization of as — 0.73 
and a n = 0.95 slope of the power spectrum. We used the 
PMTree-SPH MPI code GADGET2 (jSpringell l2005h to simu- 
late the evolution of a cosmological box with side length of 
L box = 64ft" x Mpc. 

We resimulate the region of interest in our simulation 
according to the following recipe. At z = 100, we represent 
the linear power spectrum by N ma x = 4096 3 particles of 
mass ttidm = 2.5 x 10 5 h~ x Mq. We use a Fourier transform 
of the constrained density field to substitute the overlap- 
ping Fouri er modes in to the otherwise random realization. 
Following iKlvpin et all l|200ll ). we reduce the total number 
of particles. Thus, within a lower resolution simulation of 
N = 256 3 particles, we identify the position of the model 
local group, which closely resembles that of the real Local 
Group. We then resimulate the evolution of this region of in- 
terest, using the full resolution (equivalent to 4096 3 effective 
particles) within a sphere of just 2/i _1 Mpc. Outside this re- 
gion, the simulation box is populated with lower resolution 
(i.e higher mass) particles. We are thus able to achieve high 
resolution in the region of interest, while maintaining the 
correct external environment. The initial conditions were set 
up by applying the Zeldovich approximation at z — 100 to 
avoid spurious effects due to cell crossing in the high reso- 
lution area. 

We use the same set of initial conditions to run two sim- 
ulations, one with dark matter only and another one with 
dark matter, gas dynamics, cooling, star formation and su- 
pernovae feedback. We shall call these two simualtions DM 
and SPH, respectively. In our SPH simulation, each high res- 
olution dark matter particle was replaced by an equal mass, 
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™ gas ,i "Mm e N dm 



DM: 


2.54 X 10 5 


0.15 5.29 x 10 7 


SPH : 


4.42 X 10 4 2.1 X 10 5 


0.15 5.29 x 10 7 



Table 1. Simulation parameters. From left to right, the initial 
mass per gas particle (in units of /["'Mq); mass of high resolution 
dark matter particles (in units of li"^); softening length (in 
units of h _1 kpc); Number of high resolution DM particles. 

gas - dark matter particle pair with a mass ratio of roughly 
1:5. In so doing we chose to identify the dark matter par- 
ticles in both runs with the same ID in order to simplify 
the analysis. Numerical parameters of our simulations are 
summarized in Table [1] 

For the gas dynamical SPH sim ulation, we follow the 
feedba ck and star formation rules of ISpringel fc Hernquistl 
which we briefly describe here. The interstellar 
medium (ISM) is modeled as a two phase medium com- 
posed of hot ambient gas and cold gas clouds in pressure 
equilibrium. The thermodynamic properties of the gas are 
computed in the presence of a uniform but evolving ultra- 
violet cosmic background generated from QSOs a nd AGNs 
and switched on at z — 6 (IHaardt fc Madaull 19961 ). Cooling 
rates are calculated from a mixture of a primordial plasma 
composition. No metal dependent cooling is assumed, al- 
though the gas is metal enriched due to supernovae explo- 
sions. Molecular cooling below 10 4 K is also ignored. Cold 
gas cloud formation by thermal instability, star formation, 
the evaporation of gas clouds, and the heating of ambient 
gas by supernova driven winds are assumed to all occur si- 
multaneously. 

For each star formation event, energy and metals are 
re-injected into the ISM instantaneously. Star formation is 
treated in a stochastic way and the parameters of the model 
are selected to reproduc e the Kennicut law for spiral galaxies 
l|Kenmcutdll983l , ll998t ). In order to prevent excessive star 
formation events eminating from the same SPH particle, we 
allow each gas particle to spawn at most just two generations 
of stars. Thus when a star particle is formed, the mass per 
star is half that of the progenitor gas particle, whose own 
mass is reduced accordingly. 

In order to reproduce spiral disk galaxies, we assume 
kinetic feedback in the form of strong winds driven by stel- 
lar explosions. We use the st ochastic approach developed by 
ISpringel fc Hernquistl (|2003h - SPH particles close to the star 
forming regions will participate in the wind in a probabilis- 
tic way that is proportional to the star formation rate and 
the amount of supernova energy released by massive stars. 
We assume that a fraction of /3 = 0.12 of stars will explode 
as supernovas. A fraction of this energy is used to energize 
(kick) particles in the high density regions surrounging sites 
of star formation. This removes low angular momentum gas 
from the center of dark halos, and, by inhibiting the over- 





Mft 




N p 


M b 
Mtot 


Nsats 


DM: MW : 


6.57 


230.84 


1808918 




13 


M31 : 


8.17 


248.26 


2249954 




29 


M33 : 


3.21 


181.96 


885833 




7 


SPH: MW : 
M31 : 
M33 : 


5.71 
7.81 
2.90 


220.40 
244.58 
175.71 


1745996 
2288039 
856820 


0.08 
0.12 
0.11 


20 
26 

8 



Table 2. The properties of the three local group objects, in the 
DM (top) and SPH (bottom) simulations. The colums show, from 
left to right, (1) Mtot, the total virial mass in units of 10 11 Mq, 
(2) r v ; r , the virial radius in units of kpc, (3) N p , the total number 
of DM particles within the virial radius, (4) ^ , the baryon 
fraction within the virial radius, and (5) N sa t, the number of 
satellites at z = above our 2 x 10 s Mq mass cut. 



cooling problem, contributes to the formation of extended 
gaseous disk. 

We have tested this mechanism and found it to be vital 
in the production of stable disks containing both gas and 
stars in our simulated halos. The same simuations run with 
only thermal feedback and without winds, are not able to 
create disks at all - rather the baryons and e nd up as a feau- 
turel ess spheroid of gas and stars (e.g. see lOkamoto et al.l 
120051 . and references therein). We note, that the size (ex- 
tent) of the simulated disks are related to the resolution of 
our simulations. We have run a second simulation from the 
same initial conditions and with identical star formation and 
feedback rules but with lower mass resolution (2048 3 instead 
of 4096 3 effective particles: that is, with particles that are 
8 times more massive). In this low resolution run, we find 
that the stellar disks that were produced were thicker and 
smaller. A detailed analysis of the results from these simu- 
lations will be reported in Yepes et al 2009 (in prep). 

The constrained simulations studied here are zoomed 
on a LG-like object whose three main galaxies - the Milky 
Way (MW), M31 and M33 - were reproduced with masses 
in good agreement with the latest mass estimates of these 
objects as reported in the literature (see iKlypin et al. I l2002l : 
IXueet al.ll2008l : ISeigar et al.l [20081 . and references therein). 
In Table [2] we give the characteristic properties of each local 
group halo for both our DM and SPH run. Note that the 
halos' virial radii (and therfore their virial masses) are not 
identical in the two runs. This is due to the adiabatic con- 
traction of the baryons in the SPH run which results in a 
denser core and therefore a smaller virial radius. Moreover, 
the total mass is reduced due to the lower baryon fraction 
compared with the cosmic mean. If we choose to calculate 
the mass within some fixed radius in both simulations (e.g. 
the average virial radius) and correct for the reduced baryon 
fraction, we obtain halo masses in the two simulations con- 
sistent with eachother to within a few percent. 
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2.3 Identification of Halos and Subhaloes 
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Figure 2. The average total matter density profile of our three 
halos in our DM (black) and SPH (red) simulations, normalized to 
the virial radius and virial over-density. Note that these are the 
same across much of the halo, differing only for the innermost 
regions. 



In Fig. \T\ we plot the matter distribution of the three 
main components of our local group SPH simulation. We 
plot the smoothed density of the dark matter, the gas and 
the stars. It is interesting to note that the dark matter dis- 
tribution in the DM run is the same as the dark matter 
distribution in the SPH run for the entire halo except for 
the inner most regions. Only within around 0.01 r v - u does 
the condensation of baryons in the halo's center significantly 
alter the underlying distribution of dark matter particles. 
The total matter distribution in the two runs however, dif- 
fers on a slightly larger scale. In Fig. we show the mean 
interior density profile averaged over our three halos in our 
two simulations. Across most of the halo's radius the density 
profiles are the same: only within ~ 0.1 r v i r does the total 
density profile in the SPH simulation rise significantly above 
the DM one0. 



In this section, we explain how our halo and subhalo find- 
ing algorithm works. In order to identify halos and sub- 
haloes in our simulation we have run the MPI+OpenMP 
hybrid halo finder AHF (AMIGA halo finder, to be downloaded 
freely fro m http: //popia. f t . u am. es /AMIGA) described in 
detail in iKnollmann fc Knebd (l2009h. AHF is an improve- 
ment of the MHF halo finder dGill et al1 [2004). which locates 
local over-densities in an adaptively smoothed density field 
as prospective halo centers. The local potential minima are 
computed for each of these density peaks and the gravita- 
tionally bound particles are determined. Only peaks with 
at least 20 bound particles are considered as haloes and re- 
tained for further analysis. We would like to stress that our 
halo finding algorithm automatically identifies haloes, sub- 
haloes, sub-subhaloes, etc. For more details on the mode of 
operation and actual fu nctionality we refer the read er to the 
code description paper l|Knollmann fc Knebei[2009l ). 

For each halo, we compute the virial radius r v i r , that is 
the radius r at which the density M(< r)/(4-7rr 3 /3) drops 
below AvirPback- Here pback is the cosmological background 
matter density. The threshold A v i r is computed using the 
spherical top-hat collapse model and is a function of both 
cosmological model and time. For the cosmology that we are 
using, A v i r = 355 at z = 0. 

Subhaloes are defined as haloes which lie within the 
virial radius of a more massive halo, the so-called host halo. 
As subhaloes are embedded within the density of their re- 
spective host halo, their own density profile usually shows a 
characteristic upturn at a radius r t < r v i r , where r v i r would 
be their actual (virial) radius if they were found in isolation^ 
We use this "truncation radius" rt as the outer edge of the 
subhalo and hence subhalo properties (i.e. mass, density pro- 
file, velocity dispersion, rotation curve) are calculated using 
the gravitationally bound particles inside the truncation ra- 
dius rt- For a host halo we calculate properties using the 
virial radius r v - u . 

We build merger trees by cross-correlating haloes in 
consecutive simulation outputs. For this purpose, we use 
a tool that comes with the AHF package and is called 
MergerTree. As the name suggests, it serves the purpose 
of identifying corresponding objects in the same simulation 
at different redshifts. We follow each halo (either host or 
subhalo) identified at redshift z — backwards in time, 
identifying as the main progenitor (at the previous redshift) 
the halo that both shares the most particles with the present 
halo and is closest in mass (see Eq. The latter criterion 



2 We note however, that in the analysis that follows few of our 
subhaloes penetrate the region, <~ 0.1 r v ; r , where the potential 
is significantly different in the SPH run. 



3 Please note that the actual density profile of subhaloes after the 
removal of the host's background d rops faster than for isolated 
haloes (e.g. Kazantzidis ct al. 2 0041) ; only when measured within 
the background still present will we find the characteristic upturn 
used here to define the truncation radius rt- 
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is important for subhaloes given that all their particles are 
also typically bound to the host halo, which is typically or- 
ders of magnitude more massive. Given the capabilities of 
our halo finder AHF and the appropriate construction of a 
merger tree, subhaloes will be followed correctly along their 
orbits within the environment of their respective host until 
the point where they either are tidally destroyed or directly 
merge with the host. 



2.4 Matching Halos in the two simulations 

In this section we describe how we match subhaloes in our 
SPH simulation to DM "sister" subhaloes. Our simulations 
are set up such that the identifiers (ids) of dark matter parti- 
cles are the same in both simulation. This allows us to locate 
which SPH subhalo corresponds to a specific DM subhalo. 
For a given SPH subhalo we locate the DM subhalo which 
contains both the greatest fraction of its dark matter par- 
ticles and which is closest to it in mass. We call this DM 
halo, the SPH subhalo's "sister". Specifically, for each dark 
matter halo candidate i we calculate a merit Mi, defined as 



Mi = 



(i) 



where n s h a rcd is the number of particles found in both sub- 
haloes, ridm,s is the number of particles in the DM sister 
candidate and n sp h,i is the number of particles in the SPH 
subhalo. We then identify the sister as that halo with the 
largest value of Mi . Typically, successfully matches result in 
Mi > 0.5 and we assume to have failed to find a dark matter 
sister if Mi < 0.2 for a given sph subhalo. 

We demonstrate how our sister finding algorithm works 
in Fig. |3]for one specific subhalo. The dark matter particles 
in the DM run with the same ids as the dark matter particles 
belonging to an SPH subhalo are found extended in the DM 
run. However, despite being extended, often enough of these 
particles clump together in one single subhalo, albeit with 
other dark matter particles not found in the SPH subhalo. 
We take as the DM "sister" , the entire subhalo that shares 
the most particles with the SPH subhalo. 

If we wanted to, we could have tried to find sisters by 
doing the reverse: locating particles that belong to a sub- 
structure in the DM simulation, in the SPH simulation, and 
finding which SPH subhalo hosts the most of them. The 
resulting matches are equivalent. 

It is interesting to note that our failure to match an 
SPH subhalo to a DM one is not a strong function of ra- 
dius. Although more such failures to match occur at smaller 
radii of ~ r < 0.1 r v i r ( owing to the complete de struction 
of the DM subhalo, i.e. iNagai fc Kravtsovl I2005T ). beyond 
~ 0.1 r v i r the distribution of Merit values (i.e. Eq. [TJ is uni- 
form. This implies that although the complete destruction 
of a DM subhalo may be responsible for the steepening of 
the density profile of substructure in the inner most parts 
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Figure 3. How subhaloes in the two different simulations are 
matched. On the left side we show all the dark matter particles 
belonging to a specific subhalo in the SPH simulation, in the frame 
of the host halo. The background dark matter, within which the 
subhalo is embedded, has been omitted for clarity. The red circle 
indicates the subhaloes radial extent. On the right, we show the 
same set of particles in the DM run. Note that they define an 
extended object. In red we indicate the DM subhalo to which the 
SPH subhalo was matched. Again, the background dark matter 
is omitted. Note that the DM (sister) subhalo contains additional 
dark matter particles - not found in the SPH subhalo and thus 
not plotted here - that are none the less bound to the DM sub- 
structure. 



(i.e. within r < 0.1r v i r ), it cannot explain the fact that SPH 
subhaloes are more radially concentrated across the halo's 
full virial radius. 



3 RESULTS 

In this section we explain how baryons effect the properties 
of substructures in our two simulations. Since, as mentioned 
in Sec. 12.31 we catalogue all bound objects composed of more 
than 20 particles, the vast majority of substructures in each 
of the three local group objects, have very low mass. Since 
particle mass is not constant in SPH runs (due to star forma- 
tion and feedback processes) it is incorrect to compare the 
full sub-halo population found in our DM simulation with 
that found in the SPH run. In order to compare like with 
like, we make a mass cut, comparing only the most massive 
of the z = subhaloes. We select subhaloes whose mass is 
M su b > 2 x 10 s Mq, which roughly corresponds to sub- 
haloes with more than 1000 particles. A total of around ~ 50 
subhaloes of our three main objects meet this cut in both 
simulations. This mass cut was chosen to obtain roughly the 
same number of satellites as seen orbitting about the Milky 
Way. 

In Fig. 2] we show the cumulative radial profile of sub- 
structures in our two simulations at z = 0. We make two 
different mass cuts, one considerin g the z = mass of a sub - 
structure, and a second (following Na gai fc Kravtsovll2005l ) 
which considers the subhalo's mass at infall. In order to com- 
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Figure 4. The cumulative radial density profile of substructures 
at z = 0. The thick solid and dashed line, shows this quantity 
for SPH and DM subhaloes (respectively) when considering a 
z = mass cut of M Bat ,o 

> 2 X 10 8 h _1 M Q . The thin solid 
and dashed line show this quantity for SPH and DM subhaloes 
when considering only subhaloes that were more massive then 
A^sat inf > 7 X lO 8 /i _1 M0 at infall. The infall mass cut was cho- 
sen in order to have the same number of objects as the 2 = 
cut. 



pare the same number of satellites, we consider the z — 
position of only those satellites which at infall where more 
massive then M sub > 7 x 10 s h~ I M . 

What is strikingly apparent from Fig. [4] is that the SPH 
substructures are significantly more radially concentrated 
than the DM ones. This feature is irrespective of which mass 
cut we look at, although the trend is stronger when consider- 
ing the z = mass cut. For example, within 0.5 r v i r we find 
~ 55% of SPH substructures compared with just ~ 35% DM 
ones. This result in consistent with that of (Weinberg et aD 
(2008), who compared the radial distribution of much more 
massive substructures (M su b ~ 10 10 Mq) in more massive 
systems (Mhost ~ 1O 12_14 M0) in dark matter only and SPH 
simulations of galaxy groups. 

The tendency for SPH subhaloes to be more radially 
concentrated at z — weakens with increasing redshift. Be- 
yond z = 0.75, the presence of baryons has little discernible 
effect on the radial concentration of subhaloes and the DM 
and SPH simulations return subhaloes with nearly identi- 
cal radial distributions. Indeed, as we go to higher redshifts 



Figure 5. The 2 = difference in the radial position of SPH 
satellites and their DM sisters plotted versus the radial position 
of the SPH satellite. 

we find that all SPH and DM subhaloes are less radially 
concentrated. 

Fig. [4] indicates that SPH subhaloes are more centrally 
concentrated than DM ones. Since we are able to match sub- 
haloes in both simulations, we can examine the difference in 
radial position of an individual DM - SPH subhalo pair. In 
Fig.[5]we show the difference in the radial position of two sis- 
ter halos as a function of the SPH subhalo's radial position. 
By examining how the difference in position changes with 
radius, Fig. [5] indicates that the higher radial concentration 
of SPH subhaloes seen in Fig. [4] is due to a relative deple- 
tion of DM sisters in the inner parts of the halo. Beyond 0.4 
r v i r DM sisters and SPH subhaloes are located at roughly 
the same radial distances. Within this radius however, DM 
sisters have not dissapeared (by falling below the mass cut). 
Rather, they are simply located up to 3 times further out, 
for SPH subhaloes located at ~ 0.15r v i r , and up to 1.5 times 
further out for those located at ~ 0.3r v j r . 

The scarcity of DM subhaloes in the inner parts of the 
host halo is due to the subtle changes to the substructures 
potential well caused by the presence of baryons. In par- 
ticular, the deepening of the potential well due to a cold 
baryonic component will directly impact the two main physi- 
cal processes affecting a substructure's dynamics: dynamical 
friction and tidal stripping of substructure material. 

The dynamical friction tim e scale goes as Td vn ~ Tr — » 
(e.g. lOstriker fc Tremainelll975l ; iBinnev et alj|2009h impfy- 
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Figure 6. The fraction of infall mass retained at z = for sub- 
haloes identified in SPH versus their matched DM sisters. Objects 
that sit on the diagonal line have lost the exact same amount of 
mass in both runs. Those objects in the lower triangle represent 
subhaloes which have lost more mass in the DM run while those 
objects in the upper triangular half have lost more mass in the 
SPH run. 



ing that the most massive subhaloes are those whose orbit's 
should lose the most energy and are therefore most likely to 
in spiral and be found in the inner parts of the host. Yet 
tidal stripping causes M sat to decrease, making the satellite 
less susceptible to the dynamical friction. 

The amount of material stripped can be quantified by 
looking at the ratio of the mass retained at z — to the 
mass a satellite had when it was accreted. In Fig. [6] we plot 
this ratio for SPH satellites versus the matched DM sisters. 
Rather than sitting directly on the diagonal line indicating 
equal amounts of tidal stripping, nearly all DM - SPH sub- 
halo pairs sit in the triangular lower half of Fig. [5] indicating 
that the DM sister lost more mass than the SPH subhalo 
since accretion. In some cases the difference is extreme, with 
subhaloes in SPH losing half their mass and their dark mat- 
ter sister losing 99 % of their infall mass by z = 0. Although 
a handful of SPH subhaloes have retained less than their 
dark matter sisters, the fraction of all subhaloes for which 
this is true is small (~ 5%) and the relative mass loss is not 
as extreme. 

We have seen from Fig. [6] that DM subhaloes are more 
likely to retain less of their infall mass by z = 0. Yet if our 
hypothesis is correct - that those DM subhaloes that have 



1 .0 




-0.5 i i i i i i i i i i i i i i i i i i i 
0.0 0.2 0.4 0.6 0.8 1.0 
M /M 

retained, dm retained, sph 

Figure 7. The ratio of mass retained at z = in a DM subhalo to 
that retained in an SPH halo plotted versus the 2 = difference 
in radial position of the pair. 

suffered the most tidal stripping have consequently suffered 
less dynamical friction - we should expect to see a correlation 
between the amount of mass retained and the position in 
the halo at z = 0. Indeed, those SPH halos that manage 
to fall into the centre of the halo - leaving their DM sisters 
behind in the outer regions of the halo - are also those whose 
DM sisters have lost the most mass. In other words, when a 
subhalo falls into a host halo, if it falls in early enough and if 
its orbit is such that it will lose mass due to tidal stripping, 
the stripping of the DM will decrease the dynamical friction 
on the DM subhalo, allowing the SPH subhalo to penetrate 
to the centre while leaving the DM subhalo behind. This 
notion can be tested by looking at the satellite's energy as a 
function of time. When comparing the energy loss of a SPH- 
DM satellite pair, we find that the DM subhalo has indeed 
lost more orbital energy, although the relative difference in 
energy loss is moderate (~ 10%). 

We can see this behavior explicitly in Fig. [7] where we 
plot the difference in the z — radial position between SPH 
and DM subhaloes as a function of the ratio of the infall 
mass that is retained by each subhalo z = 0. It is interesting 
to note that the :r-axis of Fig. [7] extends to ~ 1 implying 
that very few SPH subhaloes lose more mass then their DM 
sisters. The amount of mass lost by a DM subhalo is well 
correlated with the difference in z = radial position, with 
pairs that are close together having lost similar amounts of 
mass, while for those spread far apart (i.e. where ra m > 
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Figure 8. The mass of an SPH subhalo (solid line) and its DM 
sister (dashed line) as a function of redshift, normalized to the 
mass each subhalo had at infall, around z = 0.72. Note how the 
masses really start to deviate after infall, with the DM subhalo 
retaining only ~ 35% of its infall mass, while the SPH satellite 
retains ~ 55%, a factor 1.5 more. 



Figure 9. The distance from the main halo of an SPH subhalo 
(solid line) and its DM sister (dashed line) as a function of red- 
shift. Note how the trajectories start to diverge shortly after infall. 
The DM subhalo is consistently further away from the centre of 
the main halo then the SPH subhalo. 



1.5 r sp h) indicate that the DM subhalo has retained only 
< 30% of the mass retained by the SPH subhalo. 

We now follow in detail one specific SPH - DM sister 
subhalo pair, in order to highlight the salient aspects that 
are effected by tidal stripping and dynamical friction. We 
have specifically chosen a subhalo which can be used as an 
exemplary illustration for this purpose. The mass of this 
subhalo, M su b ~ 10 9 /i -1 Mq is typical of satellites in our 
sample, and was accreted at z = 0.7, an unexceptional red- 
shift. In Fig. [5]we show the mass of this subhalo as function 
of redshift (since accretion) onto the main halo. The DM sis- 
ter (dashed line) closely follows the SPH subhalo's mass for 
a short period directly after accretion. However, by z — 0.4 
the DM halo has already lost a good ~ 10% more then its 
SPH counterpart. Although the subhalo in both runs loses 
mass due to tidal stripping, by z = the DM sister has lost 
~ 65% of its infall mass, while the SPH subhalo has only 
lost ~ 45% of its infall mass. 

We may compare the mass loss presented in Fig. [8]with 
the radial distance from the centre of the main halo. In Fig.[5] 
we show the orbit of this subhalo as a function of redshift. 
Although the orbits at infall are identical, they start to di- 
verge shortly thereafter. By the first pericentric passage, the 
SPH subhalo manages to penetrate further in (by about a 



factor of 2) than the DM sister. Similarly, at apocentre the 
DM halo is ~ 0.15 r v i r further away then the SPH subhalo. 
The DM sister is consistently further behind the SPH sub- 
halo at all stages of the orbit, resulting in the SPH subhalo 
being nearly a factor of 2 closer then the DM subhalo at 

2 = 0. 

Tidal stripping is most likely to dislodge loosely 
bound material on a subhaloes periphery, where any 
change in the ho s t's po tential can have the greatest effect. 
IPefiarrubia et alj (|2008l ) showed that tidal debris fields that 
pre and succeed a subhalo on its orbit are composed mostly 
of particles depleted from this peripheral region. Since the 
cold baryons in our subhaloes have all collapsed to the cen- 
tre of the clump's potential before accretion, we except that 
the stripped material will be composed of a mixture of DM 
and hot gas particles, loosely bound to the subhalo. Indeed 
the baryon fraction of this subhalo increases as material is 
stripped from it. By z — the baryon fraction has risen by 
a factor of more than three compared with its baryon frac- 
tion at infall (not explicietely shown here). Since the total 
mass of the subhalo has decreased yet the baryon fraction 
has increased, we are led to the conclusion that the stripped 
material was mostly DM. 
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4 CONCLUSIONS AND DISCUSSION 

In this paper, we have examined the effect of dis- 
sipative baryons - treated hydrodynamically with so- 
phisticated rules for star formation and feedback - on 
halo substruc ture. We have extended the work of previ- 
ous authors (INagai fc Kravtsovl 120051 ; IMaccio et ail I2OO6I ; 
IWeinberg et aL 2008h to include not just a statistical de- 
scription of the entire subhalo population, but a detailed 
study of how a single SPH subhalo behaves differently than 
its DM counterpart. Our simulations bring a new approach 
to the study of the LG and the nearby universe. The simu- 
lations are constrained to reproduce the local environment, 
thereby enabling the testing of the processes of galaxy for- 
mation in a controlled environment that closely resembles 
the observed nearby universe. The present simulations also 
surpass previous studies in mass and force resolution. 

We use our simulations to examine the radial distri- 
bution of subhaloes. We find that subhaloes located in the 
inner parts of thier hosts and above a certain mass cut are 
more concentrated at z — in the SPH simulation than in 
the DM one. The difference in radial distribution between 
DM and SPH subhaloes above a certain mass at a given 
epoch, decreases with higher redshift, such that by z=0.75 
the two populations have similar distributions. The differ- 
ence in radial distribution weakens if subhaloes are selected 
according to infall ma ss (rather than z = mass ) but it does 
not - as suggested bv lNagai fc Kravtsovl (|2005h who looked 
at clusters 3 orders of magnitude more massive - completely 
disappear. 

The relative suppression of dark matter subhaloes lo- 
cated in the central regions of a host, is in part due 
to their susceptib i lity t o tidal forces. As suggested by 
INagai fc Kravtsovl (2005), the disintegration of DM sub- 
haloes in the inner parts of the host halo is in part respon- 
sible for their reduced numbers there. SPH subhaloes, with 
cold dense baryons in their centres, are able to inhibit the 
amount of material tidally stripped and thus hold on to more 
material. Any material that is stripped off an SPH subhalo, 
tends to originate in the outer regions of the subhalo and is 
therefore mostly composed of dark matter - little baryonic 
material is stripped. 

Since DM subhaloes lose more mass, they are also less 
likely to lose orbital energy due to dynamical friction. Our 
results suggest that although the DM subhaloes lose more 
mass due to tidal stripping, the (relatively) more massive 
SPH subhaloes spiral into the central region of the host. For 
a DM-SPH subhalo pair, this results in widely separated 
subhaloes. 

However, if tidal stripping is not particularly effective, 
dark matter subhaloes trace their SPH sisters. This princi- 
pally occurs when subhaloes are orbiting in the outer parts of 
the halo. In the inner parts, the SPH subhaloes are typically 
closer to the centre of the host than their dark matter sisters. 



The final z = masses of two sisters can differ wildly due 
to the vastly different histories each subhalo has had since 
accretion. Dark matter subhaloes almost always lose more 
mass than an SPH one, reflecting their more feeble nature. 

Our results are important for any study of substruc- 
tures in dark matter only simulations that focus on t he in - 
ner parts of the host halos. For example, IXu et all (|2009l ) 
performed an indepth study on the anomalous flux ratios 
observed with strong lensing of background QSOs. Using the 
Aquarius simulation they found that too few subhaloes exist 
in the central regions to account for the observed anomalies. 
Yet they qualify their study by explicitely stating that the 
increased survivability of subhaloes due to the presence of 
baryons will have an unknown effect on their results. Our 
work indicates that the effect of baryons may be important 
for massive subhalos that have small pericentric distances. 

Our results are important for semi-analytical models 
of galaxy formation which use merger trees extracted from 
pure dark matter sim ulations to calculate t he properties of 
satellite galaxies (e.g. lLibeskind et alJliooBT) . In these cases, 
dark matter subhaloes are often destroyed by tides earlier 
then they would be if hydrodynamics were also included. 
Thus, it is important to bear in mind that the mass of a 
DM subhalo located towards the centre, will be routinely 
underestimated. Any semi- analytical model that uses the 
mass of a subhalo extracted from a dark matter simulation 
to make predictions on the properties of the satellite galaxy 
it hosts, will suffer from this systematic effect. 

The results of INagai fc Kravtsovl <|2005h indicate that 
the bias towards more radially concentrated subhalo distri- 
butions disappears when subhaloes are selected according to 
any property that is unaffected by the processes associated 
with accretion, such as their maximum circular velocity. We 
test this hypothesis by looking at the radial distributions of 
subhaloes chosen according to infall mass (a proxy for circu- 
lar velocity) and find that although minimized (particularly 
in the inner parts of the halo where the difference in radial 
concentration is largest), the bias still exists. 

The disruption and tidal stripping of satellites has re- 
cently been addressed by a number of studies. Using con- 
trolled dark matter only simu lations and modeling the 
baryons in an idealized manner, iD'Onghia et al.l (|2009l ) ar- 
gue that the dominant mechanism in destroying substruc- 
tures in the inner regions of a galactic halo is tidal shocking. 
However, their results fail to incorporate the full adiabatic 
effect of baryons on the potential well of substructures. Be- 
ing purely dark matter, their substructures are more fragile 
and thus less able to hold onto their material. Our results 
are in agreement with Romano-Diaz et al 2008, who found 
that the number of substructures within ~ 0.1 r v i r in the 
SPH case exceeds the number in the pure DM case by almost 
a factor of two. 

Again, in the cont ext of pure dark matter simulations, 
IWetzel fc White! (|2009l ) looked at the disruption of satellites. 
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They found that the most disrupted satellites reside towards 
the centre of the halo, in agreement with this work (see. 
Fig. E| . 

Although our study is focused on the constrained simu- 
lations of the local group, it is difficult to compare our work 
with measurements of the radial concentration of the Milky 
Way's satellite population, because of the limiting biases in 
the observations. The Sloan Digital Sky Survey, which in- 
creased the number of satellites in the Mi lky Way's halo 
by a factor of two fsee IKoposov et al ] |2008l . and references 
therein) covers a just a single patch of sky. Additionally, this 
area may be biased as it points towards Virgo: a direction 
that may have an over abundance of substructures. How- 
ever, as future sky surveys such as LSST and Pan-Starrs 
reveal to us a more complete satellite populations, we will 
obtain a better understanding of the nature of the spatial 
distribution of the companions of the Milky Way. 
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